Quantum conductance of graphene nanoribbons with edge defects 
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The conductance of metallic graphene nanoribbons (GNRs) with single defects and weak disorder 
at their edges is investigated in a tight-binding model. We find that a single edge defect will in- 
duce quasi-localized states and consequently cause zero-conductance dips. The center energies and 
breadths of such dips are strongly dependent on the geometry of GNRs. Armchair GNRs are much 
more sensitive to a vacancy than zigzag GNRs, but are less sensitive to a weak scatter. More impor- 
tantly, we find that with a weak disorder, zigzag GNRs will change from metallic to semiconducting 
due to Anderson localization. But a weak disorder only slightly affects the conductance of armchair 
GNRs. The influence of edge defects on the conductance will decrease when the widths of GNRs 
increase. 
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PACS numbers: 73.63.-b, 72.10.-d, 81.05.Uw 



I. INTRODUCTION 

Recently, graphene (a single atomic layer of graphite) 
sheets were successfully isolated for the first time and 
demonstrated to be stable under ambient conditions 
by Novoselov et al.pl 13 Due to their unique two di- 
mensional(2D) honeycomb structures, their mobile elec- 
trons behave as massless Dirac fermions, 0, 0, making 
graphene an important system for fundamental physics 5, 
LJiLl' Moreover, graphene sheets have the potential to 
be sliced or lithographed to a lot of patterned graphene 
nanoribbons (GNRs) 8] to make large-scale integrated 
circuits |3 • 

The electronic property of GNRs has attracted 
increasing attention. Recent studies using tight- 
binding models [^J LUj and the Dirac equation. 12] have 
shown that GNRs can be either metallic or semi- 
conducting, depending on their shapes. This allows 
GNRs to be used as both connections and functional 
elements^]! 0] in nanodevices, which is similar to car- 
bon nanotubes(CNTs).pilG3 

However, GNRs are substantially different from CNTs 
by having two open edges at both sides (see Fig^J. These 
edges not only remove the periodic boundary condition 
along the circumference of CNTs, but also make GNRs 
more vulnerable to defects than CN TsfrU ITsl . In fact, 
nearly all observed graphene edeespl Il9ll2fj| contain lo- 
cal defects or extended disorders, while few defects are 
found in the bulk of graphene sheets. These edge de- 
fects can significantly affect the electronic properties of 
GNRs. Recent theoretical studies of perfect GNRs have 
also considered some edge corrections. ^1 |^1 |2^ But 
all edge atoms (see Fig. ^ are assumed to be identical 
and the GNRs still have translational symmetry along 
their axis. There have been no studies of changes in the 
conductance caused by local defects or extended disorder 
on the edges, which break the translational symmetry of 
GNRs. We address this issue by calculating the conduc- 
tance of metallic GNRs with such edge defects using a 
tight-binding model. 



In our calculation, external electrodes and the central 
part (sample) are assumed to be made of GNRs. And 
the edge defects are modeled by appropriate on-site (di- 
agonal) energy in the Hamiltonian of the sample. We 
utilize a quick iterative scheme [53, [24| to calculate the 
surface Green's functions of electrodes and an efficient 
recursive algorithm p3. l25j to calculate the total Green's 
function of the whole system. Finally, the conductance 
is calculated by the Landauer formula. [2rl l27j The cal- 
culation time of this method is only linearly dependent 
on the length of the sample and a GNR with disorder 
distributed over a length of 1 /jm is tractable. 
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FIG. 1: (Color online) Geometry of graphene ribbons, (a): 
A zigzag ribbon (N=4); (b): An armchair ribbon (N=7). A 
black circle denotes an edge carbon and a gray circle denotes 
a bulk carbon. A unit cell contains 2N atoms. From the top 
down, atoms in a unit cell are labeled as 1A, IB, 2A, 2B, 
NA, ATB. Atoms close to IB are 1A and 2A, and so on. 

We first study the conductance of zigzag and armchair 
GNRs with the simplest possible edge defects, a single 
vacancy or a weak scatter. Then we use a simple ID 
model to explain the zero-conductance dips caused by 
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edge defects. Finally, we study some more realistic struc- 
tures, GNRs with weak scatters randomly distributed on 
their edges. We find that a weak disorder can change 
zigzag ribbons from metallic to semiconducting, but only 
changes the conductance of armchair ribbons slightly. 
The paper is organized as follows: in section [HI we in- 
troduce the model and method employed in this paper. 
Results and discussion are presented in section ITTT1 We 
conclude our findings in section llVl 



II. 



MODEL AND METHOD 



The geometry of GNRs is shown in FigQ] A graphene 
ribbon contains two unequal sublattices, denoted by A 
and B in this paper. We use N, the number of A(B)- 
site atoms in a unit cell, to denote GNRs with different 
widths. ^3 Then the widths of ribbons with zigzag edges 
and armchair edges are W z — V3 ao and W a = -f-ao 
respectively, where ao = 2.49 A is the graphene lattice 
constant. From the top down, atoms in a unit cell will be 
labeled as 1A, IB, NA, NB. As shown in Fig. [T] an 
edge atom is a carbon atom at the edge of GNRs that is 
connected by only two other carbon atoms. In this paper, 
we consider defects that locate at the sites of these edge 
atoms only. 

The system under consideration is composed of two 
electrodes and a central part (sample). The sample (unit 
cells 1, ... , M) is a finite GNR, which may contain edge 
defects, while the left and right electrodes are assumed 
to be semi-infinite perfect GNRs. We describe the GNR 
by a tight-binding model with one 7r-electron per atom. 
The tight-binding Hamiltonian of the system is 
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where £j is the on-site energy and V pp7r is the hopping pa- 
rameter. The sum in < i,j > is restricted to the nearest- 
neighbor atoms. In the absence of defects, £j is taken 
to be zero and V pp7T =2.66 eV. [Qjj In the presence of 
defects, both the on-site energy and the hopping param- 
eter can change. Here, we only consider the variation in 
the on-site energy. A vacancy is simulated by setting its 
on-site energy to infinity (171 128|. A weak scatter caused 
by impurity or distortion will be simulated by setting 
to a small value Vt. In the case of a weak disorder, Vi is 
randomly selected from the interval ±|K. an( j om | for every 
edge atom. 

In what follows we show how to calculate the conduc- 
tance of the GNRs: 

First, the surface retarded Green's functions of the left 
and right leads (gk Q , 

l m+i) are calculated by: [23, 

il 



i = [E+I-H 0fi -Hl la A}-\ 



9o,o 
9m+i,m+i 



[E+I-H 0fi -H- lt0 A]~ 



(2) 
(3) 



where E + =E + it] (i] — > + ) 30] and / is a unit matrix. 
i?o,o is the Hamiltonian of a unit cell in the lead, and 
H— i,o is the coupling matrix between two neighbor unit 
cells in the lead. Here A and A are the appropriate trans- 
fer matrices, which can be calculated from the Hamilto- 
nian matrix elements via an iterative procedure: 0> EI3 



tnt-i to 



•tr. 



(4) 



A = t + t ti + toht 2 + ■ 
A= t + t ti + tohh + ■ ■ ■ + *o*i*2 • • • in, (5) 
where U and ti are defined via the recursion formulas 

U = (I — ti-iti-i — ti-iU-i) (6) 

U = (I — ti— i — ti-xti-x) (7) 

and 

t Q = (E + I-H 0fl )- 1 Hl lfl , (8) 

io = (E + I-H ,o)- 1 H^x,o- (9) 

The process is repeated until t n ,t n < S with S arbitrarily 
small. 

Second, including the sample as a part of the right lead 
layer by layer (from I = M to I = 2), the new surface 
Green's functions are found bv:p3l l2q 



9u = \ E+I - Hui - Hi t i +1 9i+i,i+x H ti+iY 



(10) 



Third, the total Green's function Gx,i can then be cal- 
culated by 



where 



Gxx = [E+I - Hx,x - Jy- 



^ L - ^0,130,0^0,1 

E fl = Hx,2gf^H 12 



(11) 

(12) 
(13) 



are the self energy functions due to the interaction with 
the left and right sides of the structure. From Green's 
function, the local density of states (LDOS) at site j can 
be found 
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(14) 



where Guj-\ is the matrix element of Green's function at 
site j. 

Finally, the conductance Grm of the graphene ribbon 
can be calculated using the Landauer formula[2^, 
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Here T(e) is the transmission coefficient, which can be 



(16) 
(17) 



expressed as: |32U33j 

T (E) = T I [T L Gxxr R G\ 1 ] 



where 



T L - R = i[E 



L,R 
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In this calculation, no matrix larger than N x N is 
involved. And its cost is only linearly dependent on the 
length of the GNRs. We have used this method to study 
the effects of dangling ends on the conductance of side- 
contacted CNTs.|23| We have also calculated the band 
structures of perfect GNRs by diagonalizing the Hamil- 
tonian. The conductances of perfect GNRs agree with 
the band structures. 



Zigzag ribbon, N=8 



III. 



RESULTS AND DISCUSSION 



The electronic properties of GNRs are strongly depen- 
dent on their geometry There are two basic shapes 
of regular graphene edges, namely zigzag and arm- 
chair edges, depending on the cutting direction of the 
graphene sheet (see Fig. QJ. All ribbons with zigzag 
edges (zigzag ribbons) are metallic; however, two thirds 
of ribbons with armchair edges (armchair ribbons) are 
semiconducting. [ToL Hlj . The bands of zigzag GNRs are 
partially flat around Fermi energy (Ef = eV).[lflj 
which means the group velocity of mobile electrons is 
close to zero. On the other hand, the bands of metallic 
armchair GNRs are linear around Fermi energy. 
So the group velocity of their mobile electrons around 
Fermi energy should be a constant value, which is mea- 
sured to be about 10 6 m/s.Q Since zigzag GNRs and 
armchair GNRs have so different electronic properties, 
the effects of edge defects on their conductance should 
also be very different. 



A. Single defects 

In this section, we will study the conductance and local 
density of states(LDOS) of GNRs with some single edge 
defects. A study of the effects of a single defect is not only 
realistic( e.g., a single two-atom vacancy at an armchair 
edge has been observed|l9j L but also can serve as a guide 
for us to understand the effects of more complex edge 
defects. 



1. Single vacancies 

One of the simplest defects in a zigzag GNR is a sin- 
gle vacancy caused by the loss of one edge atom (one- 
atom vacancy). In Fig|2Ja), we plot the conductance of 
a zigzag GNR (N=8) with a single 1A vacancy (dashed 
line) as a function of the energy. The solid line is for 
the perfect GNR. The defect almost does not affect the 
conductance around the Fermi energy. There are two 
conductance dips close to the first band edges and simul- 
taneously two peaks appear in the LDOS of the IB atom 
near the vacancy (dashed line in Fig|2{b)). These two 
peaks have energies different from Van Hove singularities 
of a perfect GNR, which are extreme points of the ID 
energy bands. So they are quasi-localized states caused 
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FIG. 2: (Color online) (a) Conductance of a zigzag ribbon 
(N=8, W = 17.3 A) with (dashed line) and without a vacancy 
at its edge (solid line), (b) LDOS of a IB atom when there 
is no vacancy (solid line) and when there is a vacancy on its 
nearest 1A site at the edge (dashed line). 



Armchair ribbon, N=14 
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FIG. 3: (Color online) (a) Conductance of an armchair ribbon 
(N = 14, W = 17.4 A) without vacancy (solid line), with a 
two-atom vacancy (dashed line) and with a one-atom vacancy 
(dash dot line) at its edge . (b) LDOS of a 2B atom when 
there is no vacancy (solid line) and when its nearest (1A, IB) 
pair atoms are removed (dashed line), (c) LDOS of a IB atom 
when there is no vacancy (solid line) and when its nearest 1A 
atom is removed (dashed line). 
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by the vacancy. And the conductance dips are due to the 
antiresonance of these quasi-localized states. The rela- 
tion between the conductance and quasi-localized states 
will be discussed further in Sec. IIII A 31 

The conductance and LDOS of armchair GNRs are dis- 
played in FigO In an armchair GNR, edge atoms appear 
in pairs. So a "single vacancy" can be formed by the loss 
of one edge atom (one-atom vacancy) or a pair of nearest 
edge atoms (two-atom vacancy) . These two types of va- 
cancies have very different properties especially around 
the fermi energy. For the two-atom vacancy, the conduc- 
tance is similar to that of the single vacancy situation of 
the zigzag ribbon, as well as the LDOS. However, when 
there is a one-atom vacancy, a large LDOS will be formed 
and consequently a large conductance dip will appear at 
the fermi energy. In fact, it is expected that the effect of 
a one-atom vacancy is much larger than a two-atom va- 
cancy because a one-atom vacancy breaks the symmetry 
between the two sublattices, while a two-atom vacancy 
keeps such symmetry. This is the same as in CNTs. 
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FIG. 4: (Color online) The decreasing rate of the average 
conductance (between ±0.5eV) due to a single vacancy at the 
edge. As the width of a zigzag ribbon increases, its conduc- 
tance becomes immune from a vacancy at its edge. 

The conductance around the Fermi energy is a very 
important parameter for the application of GNRs. It is 
affected by edge defects shown above. Also it depends 
on the width of the ribbon. For example, a one-atom va- 
cancy at the edge of an armchair GNR will always cause 
a zero-conductance dip at the Fermi energy. But the 
breadth of the dip will change when the width of GNR 
changes. In order to describe the effect of a defect on the 
conductance quantitatively, we introduce the decreasing 
rate of the average conductance, which is 



AG/ Go 



G(E)]dE 



J^ E G (E)dE 



(18) 



where Go(E) is the conductance of a GNR without de- 
fects, and G(E) is the conductance of the GNR with a 



defect. The decreasing rate of the average conductance 
(between ±0.5 eV) as a function of the ribbon width is 
plotted in Fig0] From the figure, we know that edge va- 
cancies affect armchair GNRs much more strongly than 
zigzag GNRs. And the effect of edge vacancies decreases 
when the width of GNRs increases. This is because there 
are more atoms in the cross section of a wider GNR, so 
electrons are easier to go around the defect. Thus we can 
use wide GNRs as connections in a nanodevice to avoid 
the change of conductance due to edge vacancies. There 
is a small bump in AG/ Go of armchair GNRs with a 
two-atom vacancy at 40A. It is because the conductance 
dips at band edges (Fig. |3J) enter into the energy range of 
±0.5 eV as the width of GNRs increases . This does not 
change the overall decreasing tendency of AG/ Go when 
the width of GNRs increases. 



2. Single weak scatters 

Another kind of single defect is a weak scatter, which 
can be caused by a local lattice distortion, an absorption 
of an impurity atom at the edge, or a substitution of a 
carbon atom by an impurity atom. Such a single weak 
scatter will be simulated by changing the on-site energy 
of an edge atom to a small defect potential V. 



Zigzag ribbon, N=8 
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FIG. 5: (Color online) (a) Conductance of a zigzag ribbon 
(N = 8) for various strengths of defect potential when the 
ribbon has a weak scatter at its edge, (b) The LDOS of an 
edge atom for V = eV (no defect) and V = 2.0 eV at that 
site. 

The conductances and LDOS of zigzag GNRs un- 
der the influences of single weak scatters with different 
strengths are presented in Fig|5] It can be seen that 
even a very weak scatter(y = 0.5 eV) can produce a 
quasi-localized state around Fermi energy and cause a 
zero-conductance dip. And the energy level and breadth 
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Armchair ribbon, N=14 
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FIG. 6: (Color online) (a) Conductance of an armchair ribbon 
(N = 14) for various strengths of defect potential when the 
ribbon has a weak scatter at its edge, (b) The LDOS of an 
edge atom for V = eV (no defect) and V = 2.0 eV at that 
site. 



of the dip increase when the defect potential increases. 
This is because the kinetic energy of mobile electrons in 
a zigzag GNR is nearly zero around Fermi energy. And 
these mobile electrons are localized to ribbon edges (edge 
states) . ^3 So they can be easily reflected by a weak scat- 
ter at the edge. On the other hand, the group velocity 
of mobile electrons in armchair GNRs around Fermi en- 
ergy is in the order of 10 6 m/s, which gives a large kinetic 
energy. So these mobile electrons will not as sensitive to 
a weak scatter as those in zigzag GNRs. In fact, there 
is no conductance dips near Fermi energy for armchair 
GNRs with a weak scatter. And the conductance dip 
at the band edge only becomes visible when the defect 
potential is larger than 2.0 eV (see FigEJ. 



3. A simple one dimensional model 

There are some common characters in the conductance 
curves and LDOS curves shown above (Figs. [3 El El 
EJ). First, there are sharp peaks in LDOS curves of per- 
fect GNRs, which are Van Hove singularities(VHS) cor- 
responding to extreme points in the energy bands. VHS 
are characteristic of the dimension of a system. In 3D 
systems, VHS are kinks due to the change in the degen- 
eracy of the available phase space, while in 2D systems, 
the VHS appear as stepwise discontinuities with increas- 
ing energy. |3J| Unique to ID systems, the VHS display 
as peaks. So GNRs are expected to exhibit sharp peaks 
in the LDOS due to the ID nature of their band struc- 
tures. Second, besides these VHS, there are new peaks 
in the LDOS of GNRs with an edge defect. And zero- 



conductance dips occur at the same energy of these new 
peaks simultaneously. These new peaks in LDOS only 
occur near the defect, but have effects on the conduc- 
tance of GNRs. So they correspond to quasi-localized 
states. And the zero-conductance dips are due to the 
anti-resonance of these quasi-localized states. The rela- 
tion between quasi-localized states and zero-conductance 
dips can be understood by a simple ID model. 
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FIG. 7: A one-dimensional model of a system including con- 
ducting bands and a quasi-localized state (QLS). The bulk of 
the system is represented by a quantum wire (QW), which 
has one conducting band. 

A GNR with an edge defect which induces a quasi- 
localized state(QLS) can be modeled as a ID quantum 
wire (QW) with a side quantum dot (see FigCJ. The 
quantum wire has one conducting band with dispersion 
relation E — 2 v cos(fcd), where E the energy of elec- 
trons, v the hopping coefficient in the QW and d the lat- 
tice spacing. The energy level of the quasi-localized state 
(side quantum dot) is El- And the coupling between the 
quasi-localized state and the QW is tx,c- If the = 0, 
the state is completely localized and has no effect on the 
conductance of the QW. When t^c ^ 0, the electrons 
not only can transport in the QW, but also can trans- 
port through "QW^QLS^QW", "QW->QLS-^QW^ 
QLS— >QW" , and so on. These different channels will 
interfere with each other and can cause resonance or an- 
tiresonance. It is easy to show that they will always cause 
antiresonance: [3^ 

To calculate the conductance of this simple system, we 
assume that the electrons are described by a plane wave 
incident from the far left with unity amplitude and a 
reflection amplitude r and at the far right by a transmis- 
sion amplitude t. So the probability amplitude to find 
the electron in the site j of the QW in the state k can be 
written as 



a k = e ikdj +re -ikdj ) ■ < Qj 



te ifc< *', j > 0. 



(19) 
(20) 



Then the transmission amplitude t and thus the conduc- 
tance of the system can be easily calculated by its tight- 
binding Hamiltonian. [35| The conductance is 
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(21) 
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From Eg 1211 we can see that when E = el, the con- 
ductance G will be zero and a dip will appear in the 
conductance curve. In other words, the incident elec- 
trons will be totally reflected when their energy is equal 
to the energy level of the quasi-localized state. So the 
quasi-localized state causes an antiresonance. This ana- 
lytical result agrees with our numerical results of GNRs. 
This relation between the conductance dips and local- 
ized states is very useful in experiments. It's not easy 
to measure the conductance of GNRs directly because of 
their small size. But the quasi-localized states can be 
find in the STS (Scanning Tunneling Spectroscopy) im- 
ages or low bias STM images, njj Then with a STS image 
or a low bias STM image, the conductance dips can be 
predicted. 



B. Weak disorders 

Experimental observed graphene edges' 0, |2(j have 
a lot of randomly distributed defects. Most of these de- 
fects are likely to be avoided in future with improvements 
in the processing of GNRs. However, as all materials have 
defects, real GNRs will always have some uncontrollable 
defects at their edges due to lattice distortion or impu- 
rity. In this section, we will consider the properties of 
GNRs under the influence of weak uniform disorders at 
their edges. An edge disorder distributed over a length 
L, will be simulated by setting the on-site energies of all 
edge atoms within a length L to energies randomly se- 
lected from the interval ±\V random \, where \V rand o m \ is 
the disorder strength. 



The conductances of zigzag GNRs with different dis- 
order strengths and distribution lengths are displayed in 
FiglHl The most important feature of the conductance 
curves is that there are gaps around the Fermi energy. 
For a N — 8 zigzag GNR with a very weak disorder 
(Vrandom = 0.25 eV) distributed over a length of 1000A, 
the conductance have a 0.25 eV gap, within which its 
maximum is less than 10 -3 • (2e 2 /h). If the disorder 
strength is 1.0 eV, the conductance gap is 1.04 eV. This is 
enormous, because a semiconducting perfect GNR with 
a similar width only has a gap less than 0.7 eV. [TlT Il2| 
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FIG. 9: (Color online) Conductance versus ribbon length 
for two zigzag ribbons (N — 8 and N = 12) with different 
disorder strengths (Vrandom = 0.25 eV and Vrandom = 1.0 
eV). The straight lines are exponential fits to the simulated 
data with the ribbon length larger than 50 A. 



Zigzag ribbon, N=8 
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FIG. 8: (Color online) Conductance versus energy for a 
zigzag ribbon (N = 8) with disorder distributed at both edges 
over a length of 100 A and 1000 A. With a weak disorder, 
zigzag GNRs change from metallic to semiconducting. 



The conductance gaps around the Fermi energy come 
from the Anderson localization of electrons. |36l I37L l38| In 
a perfect GNR or a GNR with periodic defects, the con- 
structive interference of tunneling allows that electrons 
within certain energy bands can propagate through an 
infinite GNR (Bloch tunneling). However, the disorder 
can disturb the constructive interference sufficiently to 
localize electrons. In an infinite ID system, even weak 
disorder localizes all states, yielding zero conductance. 
If the disorder is only distributed within a finite length 
L, the conductance is expected to decrease exponentially 
with length, G = Go cxp(— L/Ln), when L is much larger 
than localization length Lo-bMJ Oil We observe this is 
true in our simulations (see Fig. [5J). Each point in Fig. 
[5] is an average over several thousand disorder configu- 
rations. The localization length of electrons with energy 
close to zero is very small in the zigzag ribbons. From the 
top down, the fitted localization lengths are L = 59A, 
Lq = 53A and Lq = 45A for curves in Fig[51 respec- 
tively. So the wider the ribbon, the longer the localiza- 
tion length. And the stronger the disorder strength, the 
shorter the localization length. 

The conductance of the armchair GNRs with weak 
disorders is plotted in Fig^J Unlike the conductance 
of zigzag ribbons, there is no gap around the Fermi en- 
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Armchair ribbon, N=14 
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FIG. 10: (Color online) Conductance versus energy for an 
armchair ribbon (N — 14) with disorder distributed at its 
both edges over a length of 1000 A. 



FIG. 11: (Color online) Conductance versus ribbon length 
for two armchair ribbons (N = 8 and N = 14) with different 
disorder strengths. 



ergy. We also calculate the conductance versus disorder 
length, which is shown in Fier llll The conductance also 
decreases exponentially but much slower. The TV = 14 
armchair GNRs (W = 17.4A) have nearly the same width 
of N = 8 zigzag GNRs (W = 17.31). But their local- 
ization lengths are very different. When the disorder 
strength is V ran dom — 0.25 eV, the localization length 
of N = 14 armchair GNRs is larger than 2 /j,m, while 
L = 591 in N = 8 zigzag GNRs. So the localization 
is much weaker in armchair GNRs than in zigzag GNRs. 
As discussed in Sec. IIII Al this is because the kinetic 
energy of mobile electrons around Fermi energy in arm- 
chair GNRs are larger than in zigzag GNRs, and also 
because these mobile electrons in zigzag GNRs are local- 
ized to edges, while they distribute in the whole cross 
section of armchair GNRs. So when compared to zigzag 
GNRs, armchair GNRs are more like 2D systems where 
electrons are easier to travel around defects. There is no 
such difference between zigzag and armchair CNTs. The 
conductances of both zigzag CNTs and armchair CNTs 
are not significantly affected by disorder. [Tsj 

Recent studies of perfect GNRs have found that, with 
edge corrections which keep the translational symmetry 
of GNRs, all zigzag GNRs will be still metallic. pd| Here 
we shown that with a weak disorder at edges, zigzag 
GNRs will change from metallic to semiconducting due 
to Anderson localization. So narrow zigzag GNRs with 
a weak disorder can be used as functional elements in a 
nanodevice. This result is important because nearly all 
realistic GNRs contain some edge disorder. 



IV. CONCLUSION 

Using a tight-binding model, we have investigated 
the conductance of the zigzag and armchair graphene 
nanoribbons with single defects or weak disorder at their 



edges. We first study the simplest possible edge defects, 
a single vacancy or a weak scatter. We find that even 
these simplest defects have highly non-trivial effects. A 
single edge defect will induce quasi-localized states and 
consequently cause zero-conductance dips. And the cen- 
ter energies and breadths of such dips are strongly de- 
pendent on the geometry of GNRs. A one-atom edge 
vacancy will completely reflect electrons at Fermi energy 
in an armchair GNR, while only slightly affecting the 
transport of electrons in a zigzag GNR. The effect of 
a two-atom vacancy in an armchair GNR is similar to 
the effect of a one-atom vacancy in a zigzag GNR. A 
weak scatter can cause a quasi-localized state and con- 
sequently a zero-conductance dip near Fermi energy in a 
zigzag GNR. But its effect on the conductance of arm- 
chair ribbons near Fermi energy is negligible. The in- 
fluence of edge defects on the conductance will decrease 
when the widths of GNRs increase. Then we use a sim- 
ple one dimensional model to discuss the relation be- 
tween quasi-localized states and zero-conductance dips 
of GNRs. We find that a quasi-localized state caused by 
a defect will cause antiresonance and corresponds to a 
zero-conductance dip. 

Finally, we study some more realistic structures, GNRs 
with weak scatters randomly distributed on their edges. 
We find that with a weak disorder distributed in a finite 
length, zigzag GNRs will change from metallic to semi- 
conducting due to Anderson localization. But a weak 
disorder only slightly affects the conductance of arm- 
chair GNRs. The effect of edge disorder decreases as 
the width of GNRs increases. So narrow zigzag GNRs 
with a weak disorder can be used as functional elements 
in a nanodevice. And GNRs used as connections should 
be wider than GNRs used as functional elements. These 
results are useful for better understanding the property 
of realistic graphene nanoribbons, and will be helpful for 
designing nanodevices based on graphene. 
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